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Abstract. We use the example of the cuboctahedron, a highly frustrated molecule with 12 
sites, to study the approach to the classical limit. We compute magnetic susceptibility, specific 
heat, and magnetic cooling rate at high magnetic fields and low temperatures for different 
spin quantum numbers s. Remarkably big deviations of these quantities from their classical 
counterparts are observed even for values of s which are usually considered to be almost classical. 



1. Introduction 

Efficient low-temperature cooling by adiabatic demagnetization 
is a potential application of highly frustrated antiferromagnets 
[1-5]. However, quantitative theoretical results for the magneto- 
thermal properties have so far been obtained only in the classical 
case [1,2] and the extreme quantum case of spin s = 1/2 [3-5], 
with the exception of a study of the magnetocaloric effect in a 
(non-frustrated) ferrimagnetic spin-(5/2, 1) chain [6]. Note that 
numerical computations are already quite demanding for highly 
frustrated spin- 1/2 antiferromagnets in low dimensions and 

become quickly unfeasible for s > 1/2 or in higher dimensions. Figure 1. The cuboc- 
The cuboctahedron (see Fig. QJ is a highly frustrated tahedron: a highly frus- 
molecule which can be considered as a variant of the kagome trated magnetic molecule 
lattice with 12 sites. At high magnetic fields, it supports with 12 spin sites, 
localized magnon excitations [7, 8] which can be understood to be responsible for the enhanced 
magnetocaloric effect in these frustrated quantum spin systems. Several works have used the 
cuboctahedron as a model system for studying magneto-thermal properties [9—11]. Here we 
continue along these lines and use the cuboctahedron to investigate the behavior as a function 
of the spin quantum number s. 

2. Model, observables and method 

We consider the Heisenberg antiferromagnet in an external magnetic field h 




H = J^Si- Sj - hJ2Sf 



(1) 



where the Si are quantum spin-s operators and J > is the exchange constant between nearest 
neighbor sites The classical Heisenberg model is obtained from (JTJ) in the limit s — > oo 



at fixed J c \ = s 2 J, h c \ = sh. In this classical limit the operators Si can be considered as unit 
vectors. Accordingly, a thermodynamic quantity Q s for quantum spin s approaches the classical 
limit as 

s «o Cs (4_, M _ 2d (^ M (2) 

\S 2 J SJJ s^oo \J cl J cl J 

Using their definitions, it is easy to show that a x = = ac, cum = — 1, and ar h = 1, where %■> 
C, and M are the magnetic susceptibility, specific heat, and magnetization, respectively. Due to 
the discrete level spacing, quantum effects will be important below a crossover scale E q oc s J, 
i.e., for T < E q . At temperatures T > E q , Q c i is expected to yield an approximation to Q s 
after rescaling according to Eq. ([2]). The absolute value of the crossover scale E q is unknown, 
but one may hope that it is small in highly frustrated magnets because of the generic reduction 
of energy scales. 

We will be particularly interested in the adiabatic cooling rate at fixed entropy S 



v -I *L 
h ~T dh 



1 dM 



C dT 



(3) 



Thermodynamic quantities can be computed for the cuboctahedron by full diagonalization 
for s = 1/2 and 1. For s > 3/2 we have to resort to a truncation of the spectra to low energies 
and one should be careful about truncation artifacts. Therefore we consider only temperatures 
T < E c /10 if the spectrum is complete until an excitation energy E c even if we have partial 
results also for higher energies. For the classical case we use a combination of the Metropolis 
Monte Carlo algorithm with over-relaxation steps. 



3. Results 

Now we present some results for the low-temperature region T < s 2 J, choosing J = 1 = J c \ 
for simplicity. The behavior close to the saturation field h = 6 s is shown Figs. [2] and [3J We 
can see a trend of convergence towards the classical limit, in particular in the cooling rate Th- 
However, quantitative corrections remain large for most cases, even for spin quantum numbers 
which are as big as s = 5. The least systematic behavior at low temperatures is seen in Fig. [2] for 
h = 5 s where the system becomes gapless in the limit s — * oo. A gap opens above the saturation 
field which leads to the thermal activation of \ an d C observed in Fig. [2] for h = 7 s. Exactly 
at the saturation field, a total of 9 states with different magnetic quantum numbers becomes 
degenerate at T = [7, 10, 11]. This leads to the 1/T singularity in x an d the big values of the 
cooling rate which can be observed in Fig. [3] For s = 1/2 and h = 3 we reproduce previous 
results for the specific heat C [9]. Both the thermal activation at h = 7 s as well as the 1/T 
singularity in x at h = 6 s are quantum effects which disappear in the classical limit. 



4. Discussion 

Remarkably big deviations from the classical behavior are observed for magneto-thermal 
properties of the cuboctahedron with s < 5. This should not come as a complete surprise 
since a slow convergence towards the classical limit has already been observed for the zero-field 
magnetic susceptibility of spin chains [12] and certain magnetic molecules [13]. 

It should be noted that the cuboctahedron may not be the most representative case. Indeed, 
it has been pointed out [11] that quantum effects are more pronounced in the cuboctahedron 
than in the icosidodecahedron which can be considered as a 30-site variant of the kagome lattice. 
It would therefore be desirable to study the approach to the classical limit for non-frustrated 
lattices where Quantum Monte Carlo simulations can be used to compute magneto-thermal 
properties on bigger systems. In the meantime, one should be careful about the applicability 
of purely classical models for a quantitative description of highly frustrated antiferromagnets in 
the low-temperature region T < s 2 J. 
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Figure 2. Magneto-thermal properties slightly below (left column) and slightly above (right 
column) the saturation field: magnetic susceptibility (top), specific heat (middle), and cooling 
rate Th (bottom). All quantities are normalized to the number of spins (12). 
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Figure 3. Magneto-thermal properties 
exactly at the saturation field: magnetic 
susceptibility (top, left), specific heat (top, 
right), and cooling rate T/j (right). All 
quantities are normalized to the number of 
spins (12). 
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